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We show, for the first time, that continuous dynamical decoupling can preserve the coherence of a 
two-qubit state as it evolves during a vSWAP quantum operation. Hence, because the Heisenberg 
exchange interaction alone can be used for achieving universal quantum computation, its combina- 
tion with continuous dynamical decoupling can also make the computation robust against general 
environmental perturbations. Furthermore, since the exchange-interaction Hamiltonian is invariant 
under rotations, the same control-field arrangement used to protect a stationary quantum-memory 
state can also preserve the coherence of the driven qubits. The simplicity of the required control 
fields greatly improves prospects for an experimental realization. 
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Quantum computers use superposition and entangle- 
ment of qubits to outperform digital computers [l[ . The 
advent of these machines will unquestionably encom- 
pass a radical transformation in the way we simulate 
quantum- mechanical processes [4I, imparting a plethora 
of new achievements in science and technology. However, 
the benefits of reliable quantum information processing 
depend on the development of efficient ways to avoid or 
recover from qubit errors induced by environmental in- 
teraction Q. 

A universal set of quantum gates consists of arbitrary 
single-qubit coherent rotations and a particular entan- 
gling operation Q. Such a two-qubit unitary operation 
is the VSWAP gate, that has recently been realized ex- 
perimentally using double quantum dots [5| and neutral 
atoms in an optical lattice [y|. An ideal VSWAP gate is 
obtained by the Heisenberg coupling between two qubits, 
whose dynamics is governed by the Hamiltonian 



Ho = J<T 



(l).^(2) 
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where we use h = 1 throughout, J is the exchange con- 

( A (s) (s) (s) 

slant, and, for s = 1, 2, cr^-'-' = x.ax +yo'y +2az , where 
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ax , ay , and az arc the Pauli matrices acting on qubit 
s. Remarkably, it has been shown that the Heisenberg 
interaction alone is enough for universal quantum com- 
putation, without the need of supplementary single-qubit 
operations [7|. Thus, a protective scheme for quantum 
gates such as the vSWAP is of fundamental importance, 
since any quantum computation can be based solely on 
the exchange interaction. 

Here we show, for the first time, t he effect iveness of 
dynamical decoupling to protect the vSWAP quantum- 
gate operation on physical qubits, during a dynamical 



evolution subject to a noisy environment. Our proce- 
dure admits, not only pulsed, but also continuous con- 
trol Hamiltonians and, in fact, can culminate in a better 
protection [^. Since the Heisenberg interaction is given 
by a scalar product, it is invariant under rotations. The 
implication of such a rotation invariance is that a sim- 
ple field arrangement suffices to protect, simultaneously, 
not only the VSWAP quantum-gate operation, but also 
a quantum memory. Furthermore, the protection thus 
achieved is against general classes of errors, granting in- 
creased chances of successful experiments. 

In a previous work, we have shown that it is possible 
to use a continuously-applied external field to protect 
entangled states from errors caused by the unavoidable 
interactions between the qubit system and its environ- 
ment [9[. The question naturally arises as to whether it 
is also possible to protect an entangling operation. Here 
we show that the very same external-field configuration 
of Rcf. [9| can prevent errors from occurring during the 
application of a vSWAP quantum gate. In other words, 
we show that if the control Hamiltonian is written as 



H,(t) = n{t)- (a- 



r(l) 



r(2) 
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then the dynamics is protected by the field arrangement 
given by 



n{t) 



UzUJ [z cos (rixUJt) — y sin {n^uji)] , (3) 



where w = I-k jt^, n^ and n^ 7^ n^ are non-zero integers, 
and ic is a constant. This is a simple combination of a 
static field along the x axis and a rotating field in the yz 
plane. Moreover, addressing each qubit independently is 
not necessary; the field is supposed to be spatially uni- 
form in the neighborhood surrounding both qubits. 

The evolution operator associated with the control 
Hamiltonian is as in Ref. [9[: 
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C/,(t) = C/(2)(t)[/(i)(t) = C/(i)(t)C/(2)(f), 



(4) 



since (t^^> and a^'^' commute, where 

t/('')(i) = cxp (-icvtn^ai'A exp (-^iojin^ai^A , (5) 

for s = 1,2. Because Eq. ([1]) is a scalar product, it is 
invariant under rotations and U^{t)HoUc{t) = Hq. This 
property of the Heisenberg interaction tremendously sim- 
plifies the quantum operations executed under the pro- 
tection by continuous dynamical decoupling. If it were 
not for this rotational invariance, we would have to pro- 
ceed as in Ref. [lO| and introduce an auxiliary rotating 
reference frame, complicating the procedure. Further- 
more, this invariance has another peculiarity: the same 
field arrangement that can preserve a quantum memory, 
can also, with exactly the same configuration, protect 
the quantum-gate operation. No reconfiguration of fields 
being necessary during the gate operation is a tremen- 
dous simplification; it certainly improves the prospects 
for experimental realization. 

To illustrate our protective scheme, we begin by assum- 
ing that the interaction Hamiltonian between the qubit 
system and the rest of the universe is of the form 



iJi„, = BW.<TW+B(2) 



r(2) 



(6) 



where B*^*) 



J2m=i -Sm x„i, for s = 1,2, with Xi = X, 

(s) 

X2 = y, X3 = z, and Bm , for s = 1,2 and m = 1,2,3, 
are Hermitian operators that act on the environmental 
Hilbert space. It follows from Eqs. (|4]), ([5]), and @ that 
Uc{t) satisfies the requirement for dynamical decoupling 
El: 



UUt)HintU,{t)dt = 0, 



(7) 



where tc = 2tt/uj. In order to control the intensity of the 
exchange interaction, possible candidates for the phys- 
ical qubits should be, for example, properly built tun- 
able charge qubits [12]. Although, in this particular case, 
physical reasoning leads us to assume that each qubit is 
coupled to its own environment, we shall, for the sake 
of completeness, also study the case of a common envi- 
ronment. For our present purposes, we assume that the 
particular form of Eq. ^ is, in the case of a common 
environment, given by 



Hi, 



r(l) 



r(2; 



fAB + A*B^) 



(8) 



where B^^^ = XB + X*B\ for s = 1, 2, A is an arbitrary 
complex three-dimensional vector, and i? is a scalar oper- 
ator that acts on the environmental Hilbert space. How- 
ever, when the qubits are physically located sufficiently 
far apart, as for tunable charge qubits [T^l, it is reason- 
able to suppose that their individual surroundings act as 



uncorrelated, independent environments. In that case, 
the particular form we assume for Eq. ([5]) is written as 



Hi, 



(i).fAWi?W+A(i)*i?W 



(2).^(2)5(2)^;^(2)*^(2 



(9) 

where B^'^\ for s = 1,2, acts on the environmental 
Hilbert space of the s-th qubit, and A'*'' is an arbitrary 
complex three-dimensional vector for s = 1,2. 

We first consider the point of view of the picture de- 
fined by the unitary transformation Uc{t); thus, the total 
Hamiltonian can be written 



Hit) ^Ha + HE + Ul{t)H,,,tUc{t), 



(10) 



where He is the environmental Hamiltonian and, there- 
fore, Ul{t)HEUc{t) = He- Wc represent the environ- 
ment of each qubit as a thermal bath of harmonic oscil- 
lators. In the case of a common environment for both 
qubits, we consider He = 'Y^^uJkak^ak, where a;^ is the 
frequency of the fc-th normal mode of the common en- 
vironment, and a/c and a^^ are the annihilation and cre- 
ation operators, respectively. In the case of two indepen- 
dent and identical environments, instead of above we take 



He = ELi E/c "^kc 



(s)1'„W 



, where Uk is the frequency 



of the fc-th normal mode of the s-th qubit environment. 



and 



,(«) 



and 



,(«)^ 



are, respectively, the annihilation and 



creation operators. The frequency ujk is the same for 
both independent and identical environments. Accord- 
ingly, we take B^ = Y,k {^^dl'^^k^ + Kngka^k^'^j ,where 
gu are coupling constants. 

Starting from H(t), the Hamiltonian in the "interac- 
tion" picture is defined as 
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^^W = E E E RmAt)E^r^Kt)^(:Kt). (11) 



where a 
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^i{'{t) 



U^{t)ai'^Uo{t), for s = 1, 2 and n = 1, 2, 3, with Uo{t) = 
exp{—iHot). We have used Eq. ([6]) and defined the 
operators E^\t) = Ul{t)Bii'UE{t), for s = 1,2 and 
1,2,3, with UE{t) ~ cxp(— iiJ^t). The quanti- 






(s) 



for s = 1 , 2 and 



ties W^{t)<ji'JUc{t) 

whose matrix elements, 
Rrn,n{t), are real functions of time. We proceed as in Ref. 
[1] and assume that the absolute temperature is the same 
in the surroundings of both qubits and these qubits, as 
well as their respective environments, are identical. We 
then write down the master equation for the two-qubit 
reduced density matrix, pi(t), in the Born approxima- 
tion: 



dpijt) 
dt 



E E fdt' {vi:;::\t,t')[ai^\t),pj{t)ai^:\t')] + p^^^^^^^ ( 



s,s —1 n.n' — 1 



12) 



where we have define the coefiicients 

3 3 

■)(s.s') 



K^it^t') = E E R^nAmrn',n'{t')ci:::jit,t'), 



m—1 m' — l 

for n, n' = 1,2,3 and s, s' = 1, 2, and 



cl::2it,t') = TrE{Ei::\t)pEEi:;,\t')} 



for m,m' = 1,2,3 and s,s' = 1,2. C^;^;(i,t') is the 
correlation function between components m and m' of 
environmental operators calculated at the same qubit po- 
sition, as explained in Ref. Q. Here, Tr^ denotes the 
trace over the environmental degrees of freedom. The 
operators an (t), for s = 1,2 and n = 1,2,3, can be 
explicitly obtained as the components of the following 
vector relations: 



and 









(13) 



(14) 



where a(i) = [l + cos(4Jt)]/2, b{t) = [l-cos(4Ji)]/2, and 
c{t) = sin(4Ji)/2. The environmental density matrix. 



r 



Pe, is taken as the one for a canonical ensemble consti- 
tuting a thermal bath, that is, pE = -^ exp(—/3iJB), where 
Z is the partition function, Z = Tve [exp(— /SiJ^;)]. Here, 
/? — l/ksT, ks is the Boltzmann constant, and T is the 
absolute temperature of the environment. 
We can also write the correlation function as 

ct:Q{t^t') = T('^''^TTE{E(^\t)pEEt}{t')], 

where F*^^'* ' = 1 for the case of a single, common en- 
vironment, in which case the environmental operators 
Em (t) are independent of s, and T^'^'^ ^ = Sg^s' for the 
case of two identical, uncorrclated environments. Since 
we have 



E^nHt) - ^[A„g:4^^e— '=*-(- A:,5fc4^)^e+-'=* 



and, therefore, 



TrE{E^^Mt)PEE';:}it')} = XmX*m,Y.\9k\' 



UkC 



'iu)h(t-t') 



k 

where Uk = l/[cxp(/3wfe) — 1], we obtain 



dt 



dpiit) ^ ^ f dt'T,^-'^'\t-t')[n^^\t),pjit)[n^^'Ht')]^]+ Y^ f dt'Tt'\t-mn'^^\t)]\pi{t)n^^'\t')] 

s,s' = l-J0 s,s' = l-'0 



where 



n^'\t) 



3 3 

EE 



XmRm,n{t)^nHt): 



T}''''\t) - T(^'^">Y,\9k\'nkeM-^^kt), 



and 



Tl''''\t) 



- r(-^-^')^|gfe|2(l + nfc)exp(iwfci). 



In the limit in which the number of environmental nor- 
mal modes per unit frequency becomes infinite, we define 
a spectral density as J{uj) = J2k\9k\'^^i^ ~ ^k), with 
u! S [0,c»), and interpret the summations in 7^ (i) 
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FIG. 1: (Color Online) Amplitude Damping plus Dephasing: 
In figures (a) and (b), we show, respectively, the concurrence 
and the fidelity for independent environments and, in figures 
(c) and (d), for the common environments. The dotted (blue) 
line and solid (red) line represent the dynamics of a vSWAP 
quantum gate, with and without protection, respectively. 



and 7^^*'* '{t) as integrals over w: 
and 



, w X exp{-iojt) 
exp(puj) ~ 1 






Tl'^'\t) 



for s,s' = 1,2. Here we assume an ohmie spectral 
density with a cutoff frequency Wc, namely, J(i^) = 
?7a;exp(— cli/wc), where 77 is a dimensionless constant. To 
illustrate our method, we consider the protection of an 
entangling operation. We take J = tt/S in Eq. ^ 



and assume that the two qubits are coupled to ohmic 
environments at T = 0.2if , with the coupling constant 
f] = 1/20, and the cut-off frequency given by lOcT — 27r, 
where r = 10~^s. We consider two uncorrelated classes 
of errors: amplitude damping and dephasing. We sup- 
pose that p(0) = |tl)(tll and, in Fig. IH, we show the 
fidelities and the concurrences ^M, for the protected and 
unprotected cases of the VSWAP quantum gate. We 
consider that the qubits interact with independent or 
common environments. For the protected cases, we take 
n^ = 287r/r and n^ — IA-k/t. We observe, in all pro- 
tected cases, higher fidelities and concurrences, as com- 
pared to the unprotected cases. In all examples shown in 
Fig. ([1]), the final fidelities and concurrences of the pro- 
tected dynamics are near unity. In fact, in the protected 
cases shown, they are greater than 0.998 and higher val- 
ues can be obtained for greater values of rix and n^ . 

To summarize, we present a simplified method to pro- 
tect a vSWAP quantum gate from general classes of 
errors. Our scheme protects the logical operation at 
the same time as it is applied and, using the same 
control-field arrangement, can protect a quantum mem- 
ory or a quantum gate. The flexibility of using the 
same control field, in the static and dynamic situations, 
greatly improves the prospects for an experimental real- 
ization. Furthermore, since the quantum gates derived 
from the exchange interaction alone are universal per 
se, our methodology provides the possibility of a totally- 
protected universal quantum computation, using contin- 
uous dynamical decoupling. 

This work has been partly supported by "Fundagao de 
Amparo a Pesquisa do Estado de Sao Paulo (FAPESP)" , 
Brazil, project number 05/04105-5, and by FAPESP and 
"Conselho Nacional de Desenvolvimento Cientifico e Tec- 
nologico (CNPq)", Brazil, through the "Instituto Na- 
cional de Ciencia e Tecnologia em Informacao Quantica 
(INCT-IQ)". 



[1] D. Deutsch and R. Jozsa, Proc. R. Soc. London Ser. A 

439, 553 (1992); R W. Shor, in Proceedings, 35th An- [10] 
nual Symposium, on the Foundations of Computer Sci- 
ence, edited by S. Goldwasser (IEEE Computer Society [11] 
Press, Los Alamitos, 1994), p. 124. 

[2] R. P. Feynman, Int. J. Theor. Rhys. 21, 467 (1982). 

[3] W. H. Zurek, Phys. Today 44, 36 (1991); 
arXiv :quant-ph/0306072vl (2003). 

[4] 'd. p. DiVincenzo, Phys. Rev. A 51 1015 (1995). [12] 

[5] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 

Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and [13] 
A. C. Gossard, Science 309, 2180 (2005). 

[6] M. Anderlini, P. J. Lee, B. L. Brown, J. Sebby-Strabley, 

W. D. Phillips, and J. V. Porto, Nature 448 452 (2007). [14] 

[7] D. P. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and 
K. B. Whaley, Nature 408, 339 (2000). 

[8] J. Clausen, G. Bensky, and G. Kurizki, Phys. Rev. Lett. [15] 
104, 040401 (2010). 

[9] F. F. Fanchini and R. d. J. Napolitano, Phys. Rev. A 76, 



S. Pascazio, 
Phys. Rev. 



H. Nakazato, 
A 71, 022302 



062306 (2007). 

F. F. Fanchini, J. E. M. Hornos, and R. d. J. Napolitano, 

Phys. Rev. A 75, 022329 (2007). 

L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998); L. 

Viola, E. KniU, and S. Lloyd, Phys. Rev. Lett. 82, 2417 

(1999); P. Facchi, S. Tasaki, 

A. Tokuse, and D. A. Lidar 

(2005). 

Y. Makhlin, G. Schon, and A 

Phys. 73, 357 (2001). 

N. H. F. Shibata and Y. Takahashi, J. Stat. Phys. 17, 

171 (1977); S. Chaturvedi and J. Shibata, Z. Phys. B 35, 

297 (1979). 

M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

Cambridge, England, 2000). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 



Shnirman, Rev. Mod. 



